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SUMMARY 


Two plate elements, the linear strain triangular membrane element 
CTRIM6 and the higher order plate bending element CTRPLTl, have 
been added to NASTRAN Level 16.0. This paper describes the theo- 
retical formulation, programming details, and bulk data information 
pertaining to the addition of these elements. Sample problems 
Illustrating the use of these elements have been presented. 

A triangular shallow shell element, CTRSHL, formulated from the 
above mentioned plate elements utilizing shallow shell theory, is 
under development. Details of formulation of the element are given. 

With the addition of these elements, NASTRAN has excellent capability 
for solving plate arid shell problems using the statics, normal modes, 
and buckling analyses rigid formats. 
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INTRODUCTION 


NASTRAN presently (Level 16.0) has, in all, a total of ten different 
forms of plate elements in two different shapes (triangular and 
quadrilateral) . The present NASTRAN triangular membrane element 
TRMEM, uses a linear displacement field and is found to be overly 
stiff; the basic bending element, TRBSC, the basic unit from which 
the bending properties of the other plate elements are formed, is 
also found to be very stiff due to the constraint of normal slope 
having to vary linearly on the exterior edges of the element. The 
higher order membrane and bending plate elements (given the mnemonic 
of TRIM6 and TRPLTl, respectively) available in the literature 
(refs. 1 and 2) are found to be more accurate and computationally 
more efficient than the existing NASTRAN plate elements. In this 
paper, the addition of these elements into the Level 16 version of 
NASTRAN is described. 

A triangular shallow shell element formulated from the above mentioned 
elements utilizing shallow shell theory (mnemonic: TRSHL) is also 
being added into the Level 16.0 version of NASTRAN. Addition of this 
element is nearing completion; details of the use of this element are 
also given in this paper. 


LIST OF SYMBOLS 


a,b 

lengths 


coefficients 


constants 

D 

bending rigidity 

E 

Young ’ s modulus 

I 

moment of inertia 

L 

length 

m, n 

number of nodal lines 

M 

moment forces 

N 

membrane force 

P 

applied load 
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SYMBOLS (CONT’D) 


R 

radius 

« 

T 

temperature 

t 

thickness 

X, Y, Z 

global coordinates 

X, y, z 

element coordinates 

u, V, w 

displacements 

a, 6 

rotations 

6 

deflection 

e 

strain 

K 

rate of taper 

X 

buckling factor 

u 

Poisson's ratio 

P 

mass density 


frequency parameter 

0) 

circular frequency 

Special Symbols 

» 

denotes derivative 

Subscripts 

o 

denotes value at end of beam 

mn 

number of nodal lines 

1,2 

denotes value at top and bottom of column 

x,y 

in X and y direction 

cr 

critical 
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TRIANGULAR MEMBRANE ELEMENT, TRIM6 


The present development is based on the derivation in reference 1. 

The Important characteristics of the element are that 

1. The stresses and strains vary linearly within the element. 

2. Bilinear variation in the planar co-ordinates for the thick- 
ness of the element is permitted. 

3. Bilinear variation in the planar co-ordinates for the tem- 
perature in the element is provided. 

4. Currently the use of the element is limited to solving prob- 
lems using Rigid Formats 1 and 3. 

Element Geometry and Derivation of Matrices 

The geometry of the element is shown in figure 1. The element has 
six grid points, three at the vertices and three at the mid-points 
of the sides, u and v are components of displacements parallel to the 
X and y axes of the element co-ordinate system and are assumed to vary 
quadratically with position on the surface of the element. 


u = ai + a 2 X + a^y + a^x^ + a^xy + a^y^ ( 1 ) 

V = a? + agX + a^y + a^^x^ + a^^xy + (2) 


The undetermined coefficients a to can be related to the 12 
components of the element displacement vector. 

The thickness t of the element, and also the temperature T, at 
any point (x,y) within the element are given by 


t = 


T = 


Cl + c^x + 


d^ + d2X + 


ciy 

(3) 

diy 

(4) 
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The coefficients c- , c^, and c^, as also and d^, can be 

evaluated from the user specified 

values for t and T, respectively, at the three vertices of the 
element. 

The derivation of the element stiffness matrix, consistent mass 
matrix and consistent thermal load vector follows familiar lines 
(see, for example, ref. 3) and will not be described here. The 
stresses are evaluated at the three vertices and the centroid of 
the element. Execution time to evaluate the stiffness matrix is 

0.16 sec. per element in single precision and 0.23 sec. per element 
in double precision on the GDC 6600 machine. 

Two new bulk data cards have been added to the NASTRAN bulk data 
deck for using the TRIM6 element. They are the connection and 
property cards for the element, viz, the CTRIM6 and PTRIM6 cards, 
and are given in Appendix A. Some examples of the application of 
the element are given in Appendix B. Excellent agreement with 
plane stress elasticity and experimental solutions is shown even 
for coarse mesh subdivisions. 


TRIANGUALR PLATE ELEMENT, TRPLTl 


The mnemonic of TRPLTl will be used to refer to this element. The 
important characteristics of the element are: 

1. Bilinear variation in the planar co-ordinates for the 
thickness of the element is permitted. 

2. Bilinear variation in the planar co-ordinates for the 
temperature in the element is provided for in the 
element dependent subroutine for the calculation of 
the thermal load vector. However, this capability 
can be operational in NASTRAN only after provision 
for specification of grid point temperature gradients 
is made. 

3. Transverse shear flexibility is taken into account in 
the element stiffness formulation. 

4. Currently, the use of the element is limited to solving 
problems using Rigid Formats 1 and 3. 


Element Geometry and Derivation of Matrices 

The geometry of the element is shown in figure 2. The element has 
six grid points, three at the vertices and three at the mid-points 
of the sides. The deflection w(x,y) within the element is assumed 
to vary as a quintic polynomial in the local co-ordinates, i.e.. 


2 2 'i 

w(x,y) = + 32 + a^y + a^x + a^xy + a^y + a^x 

2 2 3 A 3 

+ 3gX y + 3gxy + a^^y + a^^x + a^^ 2 ^ y 

2 2 3 4 5 4 

+ a^^gX y + aj^^xy + a^^^y + a^^x + a^^^^ y 

4-o32^ 23^ 4 5 

+ a^gx y + Bj^^x y + a 2 Qxy + a 2 ^y'* 


( 5 ) 


The element has 18 d.o.f., viz, the transverse displacement w, 
rotation about x-axis, a, and rotation about y-axis, 6, at the six 
grid points of the element. Three constraints among the coeffi- 
cients of eq. (5) are now introduced so that the normal slope varies 
cubically along each edge. Thus the 18 independent coefficients of 
eq. (5) can be related to the 18 d.o.f. of the element. The thickness, 
t, of the element, and the temperature,!, in the element, at any point 
(x,y) within the element are given by equations (3) and (A), respec- 
tively. 

Examples of the application of the element for classical plate bend— 
ing problems are given in reference 4 and show the excellent accuracy 
of the element for such problems. Additional problems showing the 
use of the element for variable thickness plates are given in 
Appendix C. Two new bulk data cards have been added to the NASTRAN 
bulk data deck for using the TRPLTl element. They are the connection 
and property cards for the element, viz , the CTRPLTl and PTRPLTl 
cards, and are given in Appendix A. Execution time to evaluate the 
matrix is 0.987 sec. per element in single precision and 
1.579 sec. per element in double precision on the CDC 6600 machine. 

TRIANGULAR SHALLOW SHELL ELEMENT, TRSHL 

The mnemonic of TRSHL will be used to refer to this element. The 
important characteristics of the element are: 

1. Bilinear variation in the planar co-ordinates for the 
thickness of the element is permitted. 

2. Bilinear variation in the planar co-ordinates for the 
temperature in the element is provided for in the element 
dependent subroutine for calculation of the thermal load 
Vector. However, only the calculations for membrane effect 
utilize this feature at present; the calculations for the 
bending effect will be operational only after provision for 
specification of grid point temperature gradients is made. 
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Element Geometry and Derivation of Matrices 

The geometry of the element is shown in fig. 3. The shell shape is 
defined hy the height, z(x,y), above the base plane in which x,y are 
taken as local co-ordinates and X,Y,Z are the global co-ordinates. 
The shell surface over each element is approximated as a quadratic 
polynomial 



The membrane displacements u and v are assumed to vary as quadratic 
polynomials as in equations (1) and (2). The normal displacement 
w is assumed to vary as a quintic polynomial as in equation (5) . From 
the shallow shell theory of Novozhilov (ref. 5), the membrane strains 
are given by 


£ = u - z w 

xx ,x ,xx 

£ = V - Z W 

yy .y .yy 

e =u +v -2z w 

xy ,y ,x ,xy 

The element has 6 grid points and 5 d.o.f. per grid point in the 
element co-ordinate system. As in the case of TRPLTl, the three 
constraint equations are introduced so that the normal slope varies 
cubically along each edge. Additional details of formulation of 
the element and application of the element to solve shell problems 
are given in reference 6. Two new bulk data cards for the element, 
viz , CTRSHL and PTRSHL cards, that have been added to the NASTRAN 
bulk data deck are given in Appendix A. Two example problems illustra- 
ting the use of the shell element is given in Appendix D. 

CONCLUDING REMARKS 

Two plate type elements, a linear strain membrane element (TRIMS) and 
a higher order bending element (TRPLTl) , have been added to a proprie- 
tary version of NASTRAN Level 16.0. The addition of a triangular 
shallow shell element is in advanced stages of completion. The theo- 
retical formulation, bulk data information and example problems per- 
taining to these elements are discussed in this paper. The example 
problems illustrate the excellent accuracy of these elements for 
solution of plate and shell problems for statics, vibration and buckl- 
ing rigid formats. 
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APPENDIX A 


BULK DATA DECK 

Input Data Card CTRIM6 Triangular Element Connection 

Description ; Defines a linear strain triangular membrane element 
(TRIM6) of the structrual model. 


Format and Examples ; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CTRIM6 

EID 

PID 

Gl 

G2 

||E9 

G4 

G5 

G6 

+abc 

CTRIM6 

220 

666 

100 

110 

lEB 

210 

220 

320 

+C3 


s^smi 

TH 









+C3 

n 










Field Contents 


EID Element identification number (Integer > 0) 

PID Identification number of PTRIM6 property card 

(Default is EID) (Integer > 0) 


Gl, G2, G3, Grid point identification numbers of connection points 
G4, G5, G6 (Integers > 0; Gl G2 ^ G3 ^ G4 ^ G5 ^ G6) 


TH 


Material property orientation angle in degrees (Real) 
The sketch below gives the sign convention for TH. 



Remarks ; 


1. Element identification numbers must be unique with respect 
to all other element Identification numbers. 

2. Interior angles must be less than 180®. 

3. The gridpoints must be listed consecutively going around the 
perimeter in an anticlockwise direction and starting at a vertex. 


449 

























4. G2, G4, and G6 are assiimed to lie at the midpoints of the 
sides. The locations of these gridpoints (on GRID Bulk 
Data cards) are used only for global coordinate system 
definition, GPWG (weight generator module), centrifugal 
forces, and deformed structure plotting. 

5. Continuation card must be present. 
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BULK DATA DECK 


Input Data Card PTRIM6 Triangular Membrane Property 

Description ; Used to define the properties of a linear strain trian- 
gular membrane element. Referenced by the CTRIM6 card. No bending 
properties are included. 

Format and Example ; 


123456789 10 


PTRIM6 

PID 

MID 

Tl 

T3 

T5 

NSM 




PTRIM6 

666 

999 

1.17 

2.52 

3.84 

8.3 





Field 


PID 


MID 

Tl, T3, T5 
NSM 

Remarks ; 


Contents 

Property identification number (integer > 0) 

Material identification number (integer .> 0) 

Membrane thickness at the vertices of the element (Real) 
Nonstructural mass per unit area (Real) 


1. All PTRIM6 cards must have unique property identification 
nvimbers. 

2. For structural problems, the material may be MATl or MAT2. 

3. The thickness varies linearly over the triangle. If T3 or 
T5 is specified 0.0 or blank, it will be set equal to Tl. 
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BULK DATA DECK 


Input Data Card CTRPLTl Triangular Element Connection 

Description ; Defines a triangular bending element (TRPLTl) of the 
structural model* 


Format and Example ; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CTRPLTl 

EID 

PID 

Gl 

G2 

G3 

G4 

G5 

G6 

+abc 

CTRPLTl 

160 

20 

120 

10 

30 

40 

70 

110 

+ABC 


+abc 

TH 






■ 




16.2 










Field 


Contents 


FID Element identification number (Integer > 0) 

PID Identification number of PTRPLTl property card 

(Default is EID) (Integer > 0) 

Gl, G2, G3, Grid point identification numbers of connection 
G4, G5, G6 points (integer >0; Gl ^ G2 f Gh Q5 G6) 


TH Material property orientation angle in degrees (Real)- 

The sketch below gives the sign convention for TH. 


Remarks ; 

1. Element identification numbers must be unique with respect to 
all other element identification numbers. 

2. Interior angles must be less than 180®. 

3. The grid points must be listed consecutively going around the 
perimeter in an anticlockwise direction and starting at a vertex. 

4. Continuation card must be p.resent. 
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BULK DATA DECK 


Input Data Card PTRPLTl Triangular Plate Property 

Description ; Used to define the bending properties of a triangular 
plate element. Referenced by the CTRPLTl card. No membrane properties 
are Included. 


Format and Example ; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PTRPLTl 

PID 

MIDI 

11 

13 

15 

MID2 

TSl 

TS3 

4-abc 

PTRPLTl 

15 

25 

20.0 

30.0 

40. 0 

35 

3.0 

1.15 

+PQR 


+abc 

TS5 

NSM 

Zll 

Z21 

Z13 

Z23 

Z15 

Z25 


+PQR 

1.0 

9.0 

1.5 

-1.5 

2.0 

-2.0 

2.5 

-2.5 



Field 

PID 

MIDI 


Contents 

Property identification number (Integer > 0) 

Material identification number for bending (Integer > 0) 


II, 13, 15 Bending area moment of inertia pet unit width at the 
grid points Gl, G3, G5, respectively (Real > 0.0) ; 

II = Tj^^/12, 13 = 15 = Tj^/12 where T^^, T^ 

T^ are the thickness of the element at the vertices, 

respectively • 


MID2 Material identification number for transverse shear 

(Integer > 0) 

TSl, TS3, TS5 Transverse Shear thickness at the grid points Gl, G3, 
G5 , respectively . 


NSM 


Nonstructural mass per unit area (Real) 


Zll, Z21, Z13 Fiber distances for stress computation at grid points 
Z23, Z15, Z25 Gl, G3, G5, respectively, positive according to the 

right-hand sequence defined on the CTRPLTl card (Real) 

Remarks ; 


1. All PTRPLTl cards must have unique property identification 
numbers . 

2. If TSl is zero, the element is assumed to be rigid in trans- 
verse shear. 




3. If TS3 or TS5 Is 0.0 or blank, it will be set equal to TSl. 

4. IF 13 or 15 is 0.0 or blank, it will be set equal to II. 

5. The stresses at the centroid are computed at the top and bottom 
fibers. 
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? BULK DATA DECK 

Input Data Card CTRSm Triangular Shell Element Connection 

Description; Defines a triangular thin shallow shell element 
(TRSHL) of the structural model. 

Format and Example ; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CTRSHL 

EID 

PID 

Gl 

G2 

G3 

G4 

G5 

G6 

+abc 

CTRPLT 

160 

20 

120 

10 

30 

40 

70 

110 

+ABC 



TH 









S9I 

16.2 










Field Contents 


EID 


Element identification number (Integer > 0) 


PID Identification number of PTRSHL property card 

(Default is EID) (Integer > 0) 

Gl, G2, G3, Grid point identification numbers of connection 

G4, G5, G6 points (Integer >0: Gl G2 G3 G4 G5 G6) 


TH Material property orientation angle in degrees (Real) - 

The sketch below gives the sign convention for TH. 

G5 



Remarks ; 

1. Element identification numbers must be unique with respect to 
all other element identification numbers. 


2. Interior anglfes must be less than 180®. 






3. The grid points must be listed consecutively going around the 
perimeter in an anticlockwise direction and starting at a vertex. 

4. Continuation card must be present. 
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BULK DATA DECK 


Input Data Card PTRSHL Triangular Shell Property 

Description ; Used to define the bending properties of a triangular 
shell element. Referenced by the CTRSHL card. 




PID Property Identification number (Integer > 0) 

MID Material identification number for membrane effect 

(Integer > 0) 


Tl, T3, T5 Thickness for membrane action at vertices 1,3,5 of 

the elements (Real > 0.0) 

MID2 Material identification number for bending effects 

(Integer > 0) 

II, 13, 15 Area moments of inertia of the element at the vertices 

1,3,5 of the element. (Real > 0.0) 

MID3 Material identification number for transverse shear 

(Integer > 0) 

TSl, TS3, TS5 Transverse shear thickness (Real > 0.0) at the vertices 
1, 3, 5, of the element. 

NSM Non-structural mass per unit area (Real) 

Zll, Z12, Z13, Fiber distances for stress computation at grid points 

Z23, Z15, Z25 Gl, G3, G5, respectively, positive according to the 

right-hand sequence defined on the CTRSHL card (Real >0.0) 
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Remarks : 


1. All PTRSHL cards must have unique property identification numbers. 

2. If T3 or T5 equal to 0.0, or blank, they will be set equal to Tl. 

3. If 13 or 15 equal to 0.0, or blank, they will be set equal to II. 

4. If TS3 or TS5 equal to 0.0, or blank, they will be set equal to TSl. 

0.0, or blank, the element is assumed to be rigid in 
trasnverse shear. 

6. The stresses at the centroid will be computed at the top and 
bottom fibers. 
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APPENDIX B 


RESULTS OF TRIANGULAR MEffflRANE ELEMENT 
The triangular membrane element is used to solve two static beam 
problems. The first one is a cantilever beam subjected to linear 
thermal gradient over the depth as shown in figure Bl(a). The 
temperature is 188. 71-K (-120“F) at the top surface, and 322.04“K 
(120®F) at the bottom surface of the beam. This problem was ana- 
lyzed in reference 7 for demonstration of the isoparametric solid 
elements. The same material properties which were given in ref- 
erence 7 are used in the analysis. The beam was modeled with 
eight TRIM6 elements. Results for the tip displacement using 
TRIM6, IHEXi and theoretical solution are shown in table Bl. 

It is demonstrated that excellent accuracy has been obtained using 
TRIM6 with as few as eight elements. 


The second problem is a cantilever beam of variable cross section 
under a concentrated force of 133.85 lb. (595.39 N) at the free end 
as shown in figure Bl(b). The exact deflection is given by 
^y = 3.2727 ^ I ~2 § log + hx) 


- (a + bx) }] + Cj^x + I 


\ 


where 

a = 100.0 
b = 1.0 

=-(150.0 - 100 log 150) 

^2 - “(6250.0 - 100.0 (150 log 150 - 150))- 50 


(Bl) 
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Four TR1M6 elements were used to model the beam. NASTRAN and theo- 
retical results for lateral deflections along the length of the beam 
are shown In table B2. 
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TABLE Bl. TIP DEFLECTIONS FOR A CANTILEVER BEAM SUBJECTED TO LINEAR 
^ THERMAL GRADIENT :• 



Tip deflection 6 , cm 

y 

IHEXl (216 Elements) 

3.668 

IHEX2 (36 Elements) 

3.932 

IHEX3 (8 Elements) 

3.894 

TRIM6 (8 Elements) 

3.858 

Theoretical Solution 

3.762 


TABLE B2. DEFLECTIONS FOR A CANTILEVER BEAM OF VARIABLE CROSS 
SECTION UNDER A CONCENTRATED FORCE AT THE FREE END 


X 

L 

Deflection, inch 

(cm) 

NASTRAN 

TRIM6 

Exact 

0.25 

.006722 (0.0171) 

.006878 (0.0175) 

0.50 

.025601 (0.0650) 

.025658 (0.0652) 

0.75 

.053129 (0.1349) 

.053141 (0.1350) 

1.00 

.085662 (0.2176) 

.085410 (0.2169) 
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(a) Uniform cantilever beam under thermal. loading. 



0.3 IN 
(0.76 cm) 


J L 



(b) Tapered cantilever beam under concentrated load. 
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Figure B1. Beam geometry and finite element idealization 

FOR THE MEMBRANE ELEMENT (TRIM6) TEST PROBL01S. 






APPENDIX C 


RESULTS OF TRIANGULAR PLATE ELEMENT 
The triangular bending element is used to solve problems in vibration 
and buckling of thin isotropic plates. 

Vibration Problems 

I 

Two vibration problems were analyzed. The first one is a linearly 
' tapered simply supported rectangular plate shown in figure Cl (a) . 

The plate thickness is given by 

t = t^ (1 + K I ) (Cl) 

where K is a constant determining the rate of taper. Two different 
mesh .sizes of the finite element model, 1x2 and 2x4, are used. 

Due to sjnnmetry, only half of the plate is used in the analysis. 
Nondimensional fundamental frequencies for rectangular plates for three 
different aspect ratios ^ and K = 0.5 and 0.8 are presented in table Cl. 
The frequency parameter is defined as: 

o _ 2 

” (ji)d 

where co is the circular frequency, a is the length, p is the mass 
density, t^ is thickness and is the bending rigidity. Analytical 
results from reference 8 are also shown for comparison. 
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The second problem is vibration of a linearly tapered rectangular 
cantilever plate. The plate is idealized with a mesh size of 2 x 4, 
or 16 elements, as shown in figure Cl (b) . Results of frequency 
parameters as defined in equation (C2) , where m and n represent 
the number of nodal lines perpendicular and parallel to the support, 
respectively, using TRIA2 and TRPLTl are shown in table C2. Constant 
thicknesses of 0.0405 in (0.1029 cm) and 0.1215 in (0.3086 cm) was used 
when modeling with TRIA2 element. Experimental data obtained by 
Plunkett in reference 9 are also given. Tables Cl and C2 showed that 
very good results have been obtained using the higher order plate 
element. For lower modes, NASTRAN TRIA2 also gave good results. 

Buckling Problems 

The out-of-plane buckling of plate elements is evaluated from the 
differential stiffness matrix of bending plate element TRPLTl due 
to membrane prestress effects obtained from a membrane analysis 
using TRIM6 elements. To solve out-of-plane buckling of plates, a 
membrane-bending combination element is necessary. TRSHL is such a 
combination element with the added feature of membrane bending coupl- 
ing for shell problems. For pure plate problems, where the curvature 
is zero, there is no coupling between membrane and bending effects 
and for such cases TRSHL reduces to a combination element. The 
results of plate buckling in this section are obtained using TRSHL 
elements. Since the addition of TRSHL into the NASTRAN program is 
not yet complete, the results for problems in this section have been 
obtained from a stand-alone version. 
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Three buckling problems were, investigated using the triangular plate | 
and membrane elements. The Hrst one is buckling of a tapered column; 
or truncated pyramid shown in figure C2 (a). The area moment of 
inertia at any cross section can be expressed in the form 


1 = 1 . (-) 
X 1 V a / 


(C3) 


where I. is the moment of Inertia at the top of the column (x = a). 

^1 

Results for the buckling factor for a tapered column of — =0.2 

^2 

fran NASTRAN using TRIA2 and TRPLTl, and ^an analyitlcal solution from 
reference 10 are given in table C3 for comparison. Monotonic 
convergence characteristics of the buckling factor are observed 

with TRSHL element although quite slow. Results given by NASTRAN 
TRIA2 elements appear to diverge from the exact solution. The 
slow convergence of TRSHL results and the apparent divergence of 
TRIA2 results may be due to the modeling of the beam using these 
plate elements Instead of using beam type finite elements. The 
buckling factor is defined by the expression: 


El 


P = X 
cr 


(C4) 


where I 2 is the mcnaent of inertia at the lower end of the column. 

The second problem is a simply supported square plate subjected to 
uniform compression in one direction. Owing to symmetry, only one 


quarter of the plate (modeled with 2x2 mesh size) is used as shown 
in figure C2 (b) . Results of the buckling factor from NASTRAN TRIA2 
and TRSHL elements and the exact solution are shown in table C4. 

The nondimensional buckling factor X is represented by the formula 


N 

cr 


= X 



(C5) 


The third problem considered is buckling of a simply supported rec- 
tangular plate of aspect ratio -g =0.8 under inplane bending load- 
ing shown in figure C2 (c). Due to sjnmnetry, only half of the plate 
is used in the analysis. NASTRAN results using TRIA2 and TRSHL 

with different mesh sizes are shown in table C5, along with analy- 
tical results from reference 10. Table C5 clearly shows that the 
TRSHL elements gave a much better prediction of the critical buckl- 
ing load than the TRIA2 elements. 
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TABLE Cl. FUNDAMENTAL FREQUENCY FOR LINEARLY TAPERED RECTANGULAR 
PLATES SIMPLY SUPPORTED ON ALL EDGES: V = 0.3 


Aspect 

ratio 

a 

b 

NASTRAN 

TRPLTl 

Finite 

element 

layout 

Frequency parame 

2 Pt- 1/2 
iter * ua (-^ ) 

Taper rate 
K = 0.5 

Taper rate 
K = 0.8 

0.5 

1x2 

14.662 

16.242 


Theory 

15.304 

16.994 

1.0 

1x2 

24.171 

26.901 


2x4 

24.454 

— 


Theory 

24.556 

27.354 

2.0 

1x2 

58.560 

64.770 


2x4 

60.346 

— 


Theory 

60.982 

67.500 
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TABLE C2. FREQUENCY PARAMETERS FOR A LINEARLY TAPERED RECTANGULAR 
CANTILEVER PLATE; v = 0.3 


Mode 

1/2 

2 PtQ 

Frequency parameter = oj a ) 

mn mn Dq 

NASTRAN 

TRIA2 

TRPLTl 

Experiment 

m 

n 

0 

0 

2.28 

2.25 

2.47 

1 

0 

9.8 

10.0 

10.6 

0 

^B 

14.5 

13.6 

14.5 

n 


23.8 

27.0 

28.7 

B 

2 

35.9 

32.8 

34.4 

B 

3 

51.5 

47.3 

47.4 

2 

0 

31.0 

53.3 

52.5 

1 

2 

64.0 

57.7 

54.0 
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table C3. buckling factor for a tapered column 


P l^ 

Buckling factor X = -§f — 

EI2 


Finite Element layout 


TRIA2 

1.4242 

1.3618 

1.3420 

TRSHL 

1.6437 

1.6050 

1.5853 

Theory 

1.505 


TABLE C4. BUCKLING FACTOR FOR SIMPLY SUPPORTED SQUARE PLATE 
uniformly COMPRESSED IN ONE DIRECTION: V = 0.3 


Buckling factor X = 


TRIA2 

TRSHL 


Exact 


4.0356 

3.9779 


4.0000 
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1x2 Layout 


2x4 Layout 


(a) Simply supported plate. 



(b) Cantilever plate. 

Figure Cl. Plate geometry and finite element idealization for 

the TRPLTl element test problems. 





X 


(a) Taparad coluam. (b) Simply aupportad aquara plata 

undar ualf ora conpresslon. 




«o P” 


1 

1 

r t 


. 1 

1 

^ j 






(e) Simply supported plata undar Inplane banding. 


Figure C2. Column and plata geometry for TRSHL element buckling 
test problems. 
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APPEW)iX D 

RESULTS OF TRIANGULAR SHALLOW SHELL ELEMENT 

Two problems, (1) that of a spherical cap, and (2) that of a cylindri- 
cal shell roof, are considered. These are the same two example pro- 
blems analyzed in reference 11. 

Good agreement in deflections at the center of the cap is obtained 
even with relatively coarse mesh sizes as shown in table Dl. Even 
though the results appear to be oscillating about the exact value, 
the percentage error in the converged solution is very negligible. 

The geometry and the finite element idealization of the cylindrical 

shell roof used in the analysis is shown in figure Dl. Results for 

the shell roof problem and the exact solution reported by Cowper 

et. al. (ref. 11) are given in table D2. Reasonable agreement is seen 

between the finite element and the exact solutions in u., w„ and N _. 

A B xxB 

The convergence of the deflection at the central point C appears to 

be very slow; this, in turn affects the convergence of M and M 

xxC yyc* 

A non-uniform finer mesh near the center of the shell, C, is likely 

to yield converging value to the def lection w. ; the moments M ^ 

C xxC 

and Myyg may then approach the exact values. 


473 



Spherical 

Cap Problem 

Finite 

Rt/L^ = 0.02 

Rt/L^ = 0.005 

Element 

Grids 

Values of 
Etw 

c 

Values of 
Etw 

c 


pr2 

O 

P R^ 
o 

1x1 

1.151069 

1.139514 

2x2 

1.007740 

0.991776 

3x3 

1.004518 

1.001766 

4x4 

1.004365 

1.000842 


Exact 


1.009785 


1.000429 











TABLE D2. RESULTS FOR A CYLINDRICAL SHELL ROOF 
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/— N 

a *S 
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CO *ri 

*o 

0.7685 

-0.9395 

-2.0283 

-2.3828 

-2.3538 

-2.3242 

-2.0562 
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Figure Dl, GeoMtry of cylindriealehell roof and finite alaaMuit Idealization. 
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